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^_^ • Abstract 

^^ ' The Thirring model, that is, a relativistic field theory of fcrmions with a contact interaction be- 

tween vector currents, is studied for dimensionalities 2 < d < 4 using the 1/A^/ expansion, where 

r^ ■ Nf is the number of fermion species. The model is found to have no ultraviolet divergences at 

\^ • leading order provided a regularization respecting current conservation is used. Explicit 0{1/Nf) 

corrections are computed, and the model shown to be renormalizable at this order in the massless 
limit; renormalizability appears to hold to all orders due to a special case of Weinberg's theorem. 

^j. ' This implies there is a universal amplitude for four particle scattering in the asymptotic regime. 

0^ ■ Comparisons are made with both the Gross-Neveu model and QED. 
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I. Four Fermi Theories in d e (2,4) 



X 

;_( ' It has been believed for some time that a renormalizable expansion for four-fermion models exists for 



dimensions larger than two, which is naively the upper critical dimension [1 - 5]. Rather than using the 
coupling constant g^, which has inverse dimensions of mass for d > 2, to organize the expansion, the 
dimensionless parameter 1/Nf, where Nf is the number of fermion species, is used. The standard example 
is the Gross-Neveu model: 

C = i^^H^-^{i'^^^?■ (1-1) 

In this case, which has been widely studied [1-8], spontaneous fermion mass generation occurs for values of 
g^ > gc ^ 0{A'^^'^), where A is some ultraviolet cutoff in the model. It is preferable to discuss the problem 
with d a continuous parameter, d G (2, 4); the scaling properties of the model are then more transparent. Of 
course, only d — 3 can correspond to a physically realizable system. If the coupling is now fine tuned to the 
neighbourhood of gc, then light fcrmions propagate and interact via exchange of a composite scalar state of 

1 



mass 2m, where m is the fermion mass. (Actually a perfectly acceptable model also arises by approaching 
gc from, the massless phase). Because the model is strongly interacting a,t g — gc, the ultraviolet asymptotic 
behaviour of the scalar propagator, obtained by resumming a sequence of fermion - anti-fcrmion bubble 
diagrams which arc dominant at leading order in 1/iV/, is non-standard: 

^un^DAk)..^. (1.2) 

The behaviour (1.2) of Da when input to a standard power counting analysis [4,5] implies that the superfi- 
cial degree of divergence of Feynman diagrams describing corrections of higher order in 1/Nf does not depend 
on the number of interaction vertices, which in turn suggests that the expansion is exactly renormalizable. 
This property has been explicitly verified at 0{1/Nf) [5,7 - 9]. Physically, the renormalizability of the model 
may be understood as a consequence of its being the infra-red fixed point under renormalization group fiow 
of a model of fermions interacting with elementary (ie. not auxiliary) scalar fields via a Yukawa interaction 
[3,7,10]. This model is super-renormalizablc. The IR fixed point of the Yukawa model is identical to a UV 
fixed point of the Gross-Neveu model as g^ -^ g^. The relation between renormalizability and hyperscaling 
relations between the model's critical exponents (which are polynomials in 1/Nf) was stressed in [7,8]. The 
exponents are currently all known to 0{1/N'i) [11], some to 0{1/N^) [12], and have been verified for rf = 3 
by numerical simulation first for Nf = 12 [8] and most recently for Nf ~ 2 [13]. The situation is analogous 
to that in interacting scalar field theories; there the IR fixed point of the super-renormalizable 0"* theory in 
d G (2, 4) is identical to the UV fixed point of a corresponding nonlinear sigma model [1,3], which once again 
has an unexpected renormalizability in 1/Nf [14,15]. 

In this paper I wish to concentrate on another, distinct interacting fermion theory, a generalization of 
the massive Thirring model. The Lagrangian is 

£ = ^,(^+m)V;, + ;^(^,7^V*)'- (1-3) 

ZNf 

This model has also been studied in the 1/Nf expansion [16 - 19]. In this case there is no phase transi- 
tion corresponding to spontaneous mass generation, but instead the "vacuum polarization" fermion bubble 
diagrams correcting the intermediate boson propagator prove to be UV finite, despite a superficial A''^^ 
divergence, due to fermion current conservation. The situation is analogous to QED, where due to current 
conservation, the transverse projection operator {S^i^k"^ — k^k^) can always be factored from the vacuum po- 
larization, reducing the effective degree of divergence by two. In QED the result is that the photon remains 
massless at each order of perturbation theory. 

It is interesting to contrast the Thirrng model with the Gross-Neveu model. In the latter case the 
superficial A''^^ divergences do not cancel, and result in an additive renormalization of the coupling g^ , and 
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hence the need to fine-tune to recover the continuum limit. In the Thirring case, the A^^^ divergences vanish 
for essentially kinematic reasons, the coupling g^ is not renormalized (as shown below), and the continuum 
limit appears to exist, to leading order in l/Nf, for all values of g^. As in the Gross- Neveu model, the 
expression for the intermediate boson propagator in the Thirring model, which is now a vector, can be 
resummcd, to give the same asymptotic form as cqn. (1.2). Therefore the power counting arguments also 
suggest that the 1/A^/ expansion is renormalizable. However, the underlying physics is different; there is no 
phase transition corresponding to a fixed point condition, and apparently no underlying super-renormalizable 
model which yields equivalent physics in its IR limit. A proof of renormalizability, therefore, must lie entirely 
within the context of the 1/A^/ expansion of the four-fermi theory. 

The potential problem which might arise for renormalization at higher orders in four-fermi models was 
highlighted in [8]. At next-to-leading order, the corrections to the boson propagator are given by two-loop 
diagrams, exemplified in fig. 4. These are superficially A'^^^ divergent, as discussed above, but there are also 
subleading divergences of the form (on dimensional grounds) fc'^^^lnA, divergent contributions which are 
non-polynomial in external momentum, and which hence cannot be compensated by the addition of a local 
counterterm. These terms, if not removed by explicit cancellation with other divergent graphs at the same 
order (which was demonstrated for the Gross-Neveu model in [8] ) , would spoil the renormalizability of the 
model, and result in a non-local interaction being generated as the cutoff is removed. As we shall see, this 
constraint on subdivergences is a feature of a graphical expansion in which diagrams of one and two loops 
appear at the same order. It will turn out that the cancellation of non-polynomial divergences is a natural 
consequence of Weinberg's theorem [20] applied to theories with non-standard propagators. There appears 
to be no barrier to extending a renormalizability proof to all orders using standard arguments. 

The rest of the paper is organized as follows. In section II we review the work of [18,19] in setting up 
the 1/iV/ expansion for the Thirring model at leading order, extending their work, which was for the special 
(but physical!) case d = 3, to the interval d E (2,4). We shall give a closed form expression for the auxiliary 
vector propagator, and examine it in various limits, including the important deep Euclidean limit fc^ -^ oo. 
The mass of the resulting vector boson is discussed as a function of g and m. In section III the divergence 
structure of the model is discussed, 0{1/Nf) corrections computed, and the renormalization of the model 
at this order given. The condition that the fermion current is conserved translates into a requirement that 
the vacuum polarization is two-loop finite: this is verified explicitly. In section IV we compare the Thirring 
model with the Gross-Neveu model and show why the cancellation of non-polynomial divergences is to be 
expected as a result of Weinberg's theorem: the result is that both models have a very similar asymptotic 
structure corresponding to an interacting UV fixed point, despite the contrast at low energies. Finally, 



comparisons arc drawn between the Thirring model and QED, and possible implications for a non-trivial 
fixed point for the latter are discussed. 



II. Leading Order Results 

Consider the Lagrangian for the Thirring model in the following bosonised form: 

— — to — 1 

C = V'i^V'i + m-ipiipi H j==Afj,^ijfj,Tp^ + -Af^Afj,, (2.1) 

where sums on repeated spacetime indices /i and flavor indices i are understood. The field A^ is a vector 
auxiliary: it may be integrated over to recover the original Lagrangian (1.3). In d-dimensional Euclidean 
space, d G (2,4), we define {7^,7,^} = 2(5^^1, (5^^ = d, trl = 4; that is, we assume four component spinors 
and hence avoid complications due to parity violation and generation of a Chern-Simons term in d = 3 
[18,19]. The Feynman rules are thus 

fermion propagator : Sp = {i^ + m)~ 

ig (2-2) 

interaction vertex : F,, = — "/,, 

To leading order in 1/Nf the auxiliary propagator receives a contribution from "vacuum polarization", that 
is, a fermion - anti-fermion bubble (figure 1), so we write 

vector propagator : D^J(fc) = (5^^ - n^^(A:), (2.3) 

where the leading order vacuum polarization is given by 

The constant a defining the momentum routing is kept arbitrary at present. For d G (2,4) momentum 

integration is defined as follows: 

f fip') = Sd rp''-\f{p')dp, (2.5) 

Jp Jo 

with 

P^iPuIip > = -^ P f{P ) ; / PiiPvP\PK.f{p) = —^ ^.^ ^ ^^ — / p .f{p ) ; etc, (2.6) 



and 
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(4^)l"r(2 



Sd ^ ,. ..^.,, - (2.7) 



At this stage no regularization is specified. We now perform tlie trace and then apply Schwingcr parameter- 
ization: 



(2.8) 



n^,.(fc) = 4g2 / dad(3 / exp [-a{{p + akf + n?) - (3{{p + (a - l)k)^ + m^)] 
Jo Jp 

X [-2pf^p,y + 2a(l - a)kf^ky + (1 - 2a){p^^k^ +Pvk^j) + (5^^(p^ + a{a - l)fc^ + fc.p + m^)] 
Since J is now finite, the momentum p may be shifted and the integral performed; the result is 



n^.(fc) 



4g2 



dad(3 



{An) 2 Jo (a + /3)2 
2a/3 



-exp 



-{a + l3)m^ - 



2 oiP ^2 



(a + pj^ 



(a + /3) 
a/3 



d-2 



(2.9) 



(a + /3)2 2(a + /3)_ 

Note that all dependence on the momentum routing parameter a has disappeared. Now the integral of the 
second term in curly brackets, proportional to (5^^, may be reexpressed as 



dadjS 



d 



dx Jo {a + 13)2+^ x2- 



■ exp 



-X- -k — x(a + (3)m 

(a + p) 



(2.10) 



However, it can be seen that the integral in (2.10) is formally independent of x, by rescaling a and /3. Strictly, 
the integral diverges and must be made finite by use of a Pauli-Villars regulator field (eg. see [21] ch. 7). Its 
contribution to 11^,^ thus vanishes, and we can write 



n^. -Pp.(fc)n(fc2), 

with the transverse projection operator 7'^^(fc) defined as 

The remaining integrals over a and /3 are straightforward, and the result for II{k'^) is finite: 



(2.11) 



(2.12) 



n(fc^) 



,-r(2 - I) / dx- 

(47r)i ^ Jo [x{l ^ x)k^ + m^]^- 



(2.13) 



2^-fc2^(2;2-|;|; 



where F is the hypergeometric function. 

For d < 4, we see that the vacuum polarization tensor np^(fc) can be evaluated exactly at leading order, 
with the assumption of a regularization which respects current conservation. It is interesting to compare 
(2.13) with known results in d = 4 and d = 2. In the limit d ^ 4_, 

.2/1 



n^.(fc) 



67r2 V(4-d) 
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k'V^,{k), 



(2.14) 



where ^e is the Euler constant. This is almost the textbook result for one- loop vacuum polarization in 
dimensionally regularized QED4, except that since we have had no need to introduce a renormalization scale 
to make the coupling g^ dimensionless, then there is no term in InA:^. If we use the linear transformation 
properties of F to examine the limit rn? -^ (eg. [22] ch. 15), it is also possible subsequently to take the 
limit d ^ 2_(_ with the result 

n^,(fc) = -2^7'^,(fc). (2.15) 

TT 

This is exactly twice the result for the one-flavor massless Schwinger model (in that case figure 1 generates 
a dynamical photon mass g / ^/tt); the extra factor of 2 arises from our insistence on four component spinors. 

Now we set d = 3 to obtain 



o2 



n^.(fc) = -f^^M-(fc) 



1 ^(fc2-4TO2)tan-i' ^^'^' 



2(/c2)r \ 2to 



(2.15) 



This is identical to the result of Yang [18] and Gomes et al [19], with the momentum k analytically continued 
to Euclidean space. It is important to note that for the whole range d G (2, 4) the asymptotic form of Ii^,^{k) 
is not polynomial in fc^, viz: 

lim n^.(fc)^-g'Pp.(fc) ^^ i' , (2.16) 



with the numerical constant Ad given by 

_ {d-l) (4^)i 

'''^-(d-2)4r(2-f)i?(f,f-l)' ^'-''^ 

where B is the Beta function. The form (2.16) was first found by Hikami and Muta [17], modulo a difference 
of definition of trl, and a factor of d. 

Let us now return to expression (2.3) for the inverse vector propagator. Using (2.11) we can invert to 
yield the propagator 

i?..(fc)=7',.(fc)^-^ + ^. (2.18) 

As argued in [16,18,19], the second term in I?^^(/c), which is longitudinal, behaves as a constant in the limit 
k — > 00, and might naively be expected to lead to poor ultraviolet behaviour. However, since the vector 
auxiliary interacts with a physical current which is conserved, S-matrix elements and observables constructed 
as gauge-invariant combinations of the fields "0 and A^ (in the sense used in QED) should not be affected by 
this problem, although Green functions in general might be. We shall see this when we calculate 0{1/Nf) 
corrections in the next section. To aid calculation (but not to define the vector propagator, as would be the 
case in QED), following [19], we introduce a gauge-fixing term -^{df^A^)"^ to the Lagrangian (2.1), which 

6 



has the effect of moderating the UV behaviour of the second term for finite fi, which has dimensions of mass. 
The vector propagator becomes 

^.^(^;m)-^.^Wt3^ + ;;^^- (2-19) 

The scale fj, is in effect a regulator which should not appear in final expressions. The limit fi -^ oo recovers 
the original Thirring model, whereas the limit ji ^ specifies a "Landau gauge" which will render the 
two-loop calculation in the next section much easier. Note also that in the limit fc^ -^ 0, Il(k'^) -^ and 
Dfii,{k) -^ S^ij: hence the infrared problems associated with QED are not present here. 

On the assumption that the longitudinal piece of Z?^^ has no physical consequence, we focus on the 
transverse piece and identify a pole condition for the mass of the vector My : 

l-g^Ml ^\ ^^ F(2:2-|;|:P^)^0. (2.20) 

In general this is a transcendental equation. It can be solved in two limits. For strong coupling g^ 3> ra^^'^ the 
vector channel will be dominated by a tightly bound fermion - anti-fermion state, so My <C m?- Therefore 
we can expand F to obtain 



A4 _ m^ ( 3(47r)i 
^2~ TFTo i^ 



(2.21) 



m^ r y4r(2- 2; / 

For arbitrarily weak coupling, real solutions of (2.20) can only exist if the hypergeometric function is able 
to grow arbitrarily large. We expect a weakly bound state to have mass given by 

Mv = 2m- £. (2.22) 

As Mv -^ (2?7i,)_, the hypergeometric function diverges only for d < 3. In this case we can once again 
perform a linear transformation on F to get 

l-.^M^^^^f^l'^fl + oO-O, (2.23) 



TO4-d(47r)2 Vm 



m 



ie. the binding energy is 



2^d-2 ^\ 2 > 



2 
3-d \ \ 3-d 



e = m\g'm''-' \\U ■ (2-24) 

The case d = 3 must be handled separately; we use expression (2.15) to find the binding energy is essentially 
singular in g^: 

e = 2mexp( ^). (2.25) 
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Finally, for d > 3 the hypergeoinctric function remains finite as fc^ -^ —Am'^. In this case the bound state 
vanishes (ie. the would-be pole coalesces with the branch cut in F) for values of g below a critical gc given 

by 

,_, (4^)i(rf-l)(d-3) 

dr — ^ : t: • (2.2b) 

^^ 16r(2-|) ^ ' 

In the subcritical region Df^i, has no poles on the physical sheet; the vector can only be regarded as a resonant 

intermediate state in four fermion scattering. 

In the deep Euclidean region k'^ ^ oo things simplify considerably: the vector propagator has the form 

hm D^,{k)=V^4k)^—^. (2.27) 

In this limit the four fermion scattering amplitude has the form AdJfi{q)Jfi{q + k)/Nfk'^^'^. As we shall argue 
in the final section, this interaction receives no corrections in the l/Nf expansion, and is thus a universal 
form characterizing the short distance structure of the model; in other words it defines a UV fixed point. 
In this respect it resembles the Gross-Neveu model as discussed in [8] (though note the definition of A^ is 
distinct). In the next section when the rcnormalization of the model at 0{l/Nf) is discussed, the form (2.27) 
will be used throughout. 



III. Renormalization at 0{l/Nf) 

In this section I will discuss the renormalization of the model to next-to-leading order in the 1/A^/ 
expansion. First let me review why we might expect such a programme to be feasible. The short distance 
fluctuations of the model are encoded in the asymptotic form for the vector propagator (2.27). Suppose we 
analyze the superficial degree of divergence of a higher order diagram with N^ external fermion lines, Na 
external auxiliary lines, and V vertices. If we use the form (2.27), then standard power counting analysis 
gives the superficial degree of divergence uj: 

uj^d-^-j^N^-NA. (3.1) 

It is interesting to compare this result with that for canonical boson asymptotics, D{k) ^ 1/fc^, which 

applies, say, for QED: 

d-1 d-2 4-d , , 

tocan ^d —N^ —Na —V. (3.2) 

For d < 4 the degree of divergence falls as the number of vertices increases: this is characteristic of a 
super-renormalizable theory. Only when uj is independent of V can a perturbative expansion be exactly 
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renormalizable (ie. divergent graphs appear at every order of the expansion, but can always be made finite 
by retuning a finite set of counterternis). 

Using (3.f ) we can compile a list of potentially dangerous graphs for d E (2, 4). The 0{l/Nf) contributions 
are shown in figures 2-4. The fermion self-energy ip-tp (fig. 2) has lu ~ 1, as in QED4, but since that 
divergence is odd in loop momentum, the true divergence is logarithmic {lo = 0). The vertex correction ipipA 
(fig. 3) also has uj = 0, but the four-vector AAAA scattering, which is superficially divergent in QED4, here 
has UJ ^ d — 4 and so is safe. One-point and three-point vector scattering vanish by Furry's theorem, leaving 
the vector two-point function AA (fig. 4), with lo = d — 2 as the only other superficially divergent case. One 
can then argue [19] that in a regularization which respects current conservation, one can always extract a 
factor k'^Vfiv{k) from these diagrams to give uj = d — 4, and hence no new divergence. The model (2.1) can 
then be renormalized simply by rescaling the tp and A fields, and retuning the fermion mass. I shall show 
in this section that this conclusion is correct, though the argument is not quite so straightforward. Because 
the vacuum polarization is non-polynomial in fc^, it is in fact only permissible to extract 'P^i,i,(k) from inside 
the graph, which docs not improve the power-counting. As we shall see, there are divergent contributions 
both of degree lo = d—2 and lj ^ Q. The applicability of power counting to four fermion models is discussed 
further in the final section. 

It is worth contrasting the four-fermi case with the situation in pure scalar theories. In the renormalization 
of the nonlinear sigma model in d e (2,4) [15], the power counting gives 

LO = d- -^N^ " 2Na, (3.3) 

where now tp denotes the elementary scalar field and A the auxiliary scalar boson. The only superficial diver- 
gences are ipip {lo — 2) and ip-ipA (a> = 0); the auxiliary propagator AA has lo = d— 4 and hence is superficially 
convergent. On dimensional grounds there is no reason to expect any non-polynomial divergences. 

The procedure for renormalizing the model follows the treatment in [8]: first we redefine the Lagrangian 

£ - Z^^,(^ + M)^, + ^|=Z^z|^,4 V^ + ]:Za ( AI + -^{d^A^^f] . (3.4) 

The constants Z-^, Za, M and in principle g and fi are all cutoff-dependent, and must be adjusted at each 



order of the l/Nf expansion to keep physical matrix elements finite. As we have seen, at leading order an 
adequate choice is: Z^ — 1; M ^ m, the physical fermion mass; fi -^ 00; and g, Za unconstrained. Za 
simply defines the scale of an auxiliary field at leading order and hence has no physical relevance. The first 
divergent Green function we must examine is the fermion self-energy (fig. 2): 
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Note that with the definition (3.4), the vector propagator at leading order is Z j^ Df^i,{k; ^), with £)^^ defined 
by (2.19). On rearranging, we find 






A{p^)V,Ap) + B{p';f.)^ 



with 



hm A(A:2) = ^1^— 



; B{k^) 



(3.6) 



(3.7) 



We will treat the parts depending on A(p'^) and B{p^] p) separately. For the first piece, apart from a term 
which is odd in p and hence vanishes on J , the leading contribution is 0{p^^dp), and hence logarithmically 
divergent. With the choice of a simple momentum cutoff \p\ < A, and the definitions (2.5,6), we find 



S^(fc) = -Z, 



Cd {d - If 

Nf2{d^2) 



(d-4) 



M 



In — + finite. 



with the constant C'd defined, as in [8]: 



C,= 



B(f,2-|)i?(|,|-l) 



(3.8) 



(3.9) 



Note Cd is positive definite for d E (2, 4). For d = 3, Cd has the value 4/7r^. The piece depending on B(p^; p) 



can be recast in the form 

j.i3(^) ^ _ .g^y r(2~|) /•! ^^^^^^ _ ^^^2 ^ ^^2 ^ (^ _ a;)/i2]i-2 



^/ (4^)^ Jo 



2i 

i(l + x)^ +M- — 1^ 
a 



In the limit p^ ^ oo: 



j.B(^) ^ _ g'^v./ 'r(2 J) 2 ^^^ ^ ^^^ ^ o(m''-*). 



+ o(m''"')- 

(3.10) 
(3.11) 



Nf (4^)i (d-2) 

Note that S'^ and S^ have different characteristics: E"^ is cutoff-dependent but independent of the coupling 
constant 17, whereas S^ is finite for finite /i, but depends on both p and g. The Thirring model limit /i^ ^ 00 
cannot be taken for S in isolation. 

We can now write for the full inverse fermion propagator 

Sp'^{k) = Z^{ilf: + M - S(fc)) = i^ + m, (3.12) 

thus defining the wavefunction renormalization constant Z^ and the renormalized (ie. physical) mass m in 
terms of the bare mass M: 



Z^ — 1 



Cd{d-lf{d-^)^^K .gV-' 2r(2-|) 



In — 



m^ M [1 



Nf 2d(d-2) m Nf (47r)i(d-2) 

Cd2{d-l)\ A 



finite. 



Nfdid-2)'^^+''^''V^''- 



(3.13) 
(3.14) 
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Thus wc find an expression for tlie pliysical mass m independent of the regulator fi, and a wavefunction 
constant Z^ which depends on /i. This is, of course, very similar to what is found for QED: the physical 
parameter m is renormalized in a gauge-invariant fashion, whereas the unphysical Z-^ is not. The term 
proportional to C'd in (3.13) was derived in [17], but the /i-dependent piece was neglected. 

Next we calculate the 0(1 /Nf) contribution Pj^ to the vertex (fig. 3). For zero external momentum we 

have 

^[1] ig^ Z^Zl r 1 1 n , X 

j^2 Z^Za Jp '^ip+M ip+M 

i 
Z I Z^ a^ f 1 



^ Nf Jp{p^ + M^)^^ 

where the second line follows because we are only interested in the divergent part. Using the same procedure 
as before, we find for the full vertex 



A 


— ^ X 


+ rw 






= -ig 


1 
z^zx 


ix 



— Cd- 7-^ — ^ — - In h g u -, 2-^ + finite 

Nf I 2d{d-2) m ^ ^ {A7r)i{d-2) 



(3.16) 



However, the constant Z^ has already been determined in (3.13), and is found to exactly cancel both A- and 
/i-dependent terms in (3.16). Hence 

Ta = r=ZXlx [1 + 0{1/Nf) X finite] . (3.17) 

Once again, this is a familiar situation from QED: current conservation plus gauge invariance ensures that 
the divergent and gauge-dependent parts of the self-energy and vertex corrections cancel, ie. Zi = Zi (as 
noted in section H, there are no problems with infrared divergences in the Thirring case). We expect this 
cancellation to persist at higher orders. So, to maintain the finiteness of fermion self-energy and vertex 
corrections to 0{1/Nf), our only requirement of Za and g so far is that the combination Z^g be finite. 
However, we have not yet exhausted the list of superficially divergent graphs. We next consider the 0{1/Nf) 
corrections to the vector propagator, which consists of two two-loop diagrams (figs. 4a,b): 

nW (fc) - -^ X |||^ X Nf X [2/^-,(fc) + Ilik)] ; 

..,._/■/ 1 1 1 1 

The fermion masses, which have little impact on the ultraviolet behaviour of the integrands, have been 
neglected (see below). Now, given the asymptotic form (2.27) for Dap(q) ^ Ad / g'^ q'^^'^ , it is easy to see 
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that the combination of rcnornialization constants multiplying /";/ reduces to ZaQ^, which from previous 
considerations must be cutoff-independent. The conclusion is that n[,i, must be UV finite if the model is 
to be consistently renormalized at this order. The constant Za, being just the scale of an auxiliary field, is 
undetermined in the model, and the value of g also appears to be irrelevant as regards the UV behaviour of 
the model (Yang [18] points out that g cannot be renormalized, since it appears in the Lagrangian (3.4) in 
a gauge- variant manner after rescaling A^ h^ Afj_/g). 

Now, as shown in [8], diagrams such as those of fig. 4 generically have divergences of two forms, one 
independent of fc, proportional to A''^^ + const x ikf'^^lnA, and the other proportional to fc'^^^lnA. We 
must show that for the Thirring model both occur with zero coefficient. First we deal with the momentum- 
independent piece, following the technique used in the appendix of [5]. 

Consider the expression (2.4) for the one- loop vacuum polarization Iv^l{k), and in particular the result 
of differentiating it with respect to external momentum k. By using the identity 

1 I = I . I (3 19) 

we see that each differentiation is equivalent to a zero momentum insertion of —iA^ (modulo a factor of 
g/^/N}). Thus 

gl^n™ (fc) = -(2(1 - af + 2a2) j;,„^(fc) + 2a(l - a) J^,„^(fc), (3.20) 

where J^^l^Bi^) ^^^ represented in fig. 5. However, as we showed in section II, IlJ^L and hence J"j^„«, are 
finite analytic functions of k which are independent of the momentum routing a. Hence 

^ ^al^ng(fc) = 2(1 - 2a)[2J,^„,(fc) + Jl^pik)]. (3.21) 

Now we contract the R.H.S. of (3.21) with Dap{k) and perform J^, to obtain the two-loop integrals /"^j* at 
zero momentum. Since the R.H.S. must vanish for arbitrary a, we conclude 

2/^.(0) +/^(0) = 0. (3.22) 

This argument shows that the momentum-independent part oiIi\il{k) is identically zero, and hence that we 
do not need to worry about A''^^ divergences. However, each diagram is individually divergent, as can be 
seen by setting a = in (3.20), then evaluating J^, Jai,api^)-^ai3{k); we find for the leading divergence 

21V(0)-V(4^^.^(,^ d(rf-2) ' ^^-'^^ 

where A is a momentum cutoff for J, . 
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Of course, we might have anticipated that the A^^^ superficial divergence vanishes due to current conser- 
vation, as it did at leading order, and indeed as it does in QED to all orders, where 11^^(0) = ensures that 
the photon propagator retains a zero mass pole, as required by gauge invariance. However, no such argument 
constrains the momentum-dependent divergence fc'^^^lnA - for instance, in QED4, /c^lnA divergences are 
physical and responsible for charge renormalization - and it is not a priori clear what will happen in a model 
with non-standard asymptotics in d < 4. To examine this case it is necessary to perform an explicit two-loop 
calculation, using the techniques developed in [8], which are now outlined. 

After performing the trace over Dirac indices, each integral may be reduced to several components of the 
form 

constantx/// ^7^^^'^^ 1"^°^^ "\ (3.24) 

using Feynman parameterization and momentum shift in p. Here P and Q are polynomial functions of 
momenta, Xi are Feynman parameters which are to be integrated over a domain S^x^ < 1, and n is integer. 
The algebra is considerably simplified by the choice of "Landau gauge" p^ -^ 0, which means that all terms 
proportional to q^, qp can be discarded. The integral over the fermion loop momentum p can always be 
performed for n > 2: 



lp[p^ + QY r(n)(47r)i 

Momentum dependent divergences arise when there are two or more Feynman parameters, in which case 
there is an intermediate integral of the form 



i-y.,^.x, ^ A + Bx 

(a + bxi + cxj) 



<^ r,r:\::J^r ^ (3-26) 



where s is non-integer, and the coefficients 6, c are O(g^), where q is the remaining loop momentum, but the 
coefficient a is O(q'^k^). The contribution to (3.26) from the Xi ^ limit of the integral is then 

Ab ( (3 - 2s) 2ac\ B ^ / 1 \ , , 



(s-l)Aa''-iV (2-s) A J (s-l)(s-2)Aa''-2 ^^6^ 

with A = 4ac — b^. The exact expression on which this approximation is based is given in the appendix 
of [8]. Note that all factors, including A, must be expanded consistently to 0{l/q'^) due to the presence of 
0{q^) terms in the numerator polynomial P. Once this limit has been isolated, the remaining integral over 
q is of the form 

| i?.,..(M)fc^-^ (3.28) 

Jq r 

where R is 0{k'^q^). The form (2.27) for Da(3 is now sufficient to evaluate J with a momentum cutoff: 
it yields a logarithmic divergence. Any remaining integrals over Feynman parameters give combinations of 
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Beta functions in d. The final result is 



2n[H(fc) = -n[V^i(fc) 



~Wf {Anrnd) d ^A.-(^)(^ )' l^X (3.29) 

2^^ n,.ik^)^~WCdid-l)Hd-4)^^A 



on rearranging. We note that each diagram is individually transverse and equal and opposite respectively 
to twice the /i -independent parts of the self-energy correction (3.13) (a), or the vertex correction (3.16) (6). 
The main result, of course, is that the two- loop vacuum polarization Ir^l{k) is UV finite in the massless 
theory. Strictly, we have not demonstrated the independence of this result on the "gauge fixing" parameter 
fi, and must rely on reasoning that the vacuum polarization, which yields charge renormalization in QED, 
is gauge invariant. Perhaps more importantly, we have not considered divergences of the form Af^^^fc^ In A: 
these probably exist and correct the mass of the vector My following the discussion (2.20-26). 

We have now exhausted the list of divergent Green functions at 0{l/Nf), and proven that the model is 
renormalizable at this order (at least in the massless limit) - indeed the only physical (ie. gauge invariant) 
renormalization that must be made is that of the fermion mass (3.14). 



IV. Discussion 

We begin the final section by thinking about why the cancellation of k'^^'^ In A divergences in the two-loop 
vacuum polarization takes place. The way the calculation has been presented here suggests a cancellation 
between the two diagrams n[i^ and Uliu (fig. 4). However, a similar analysis of the Gross-Neveu model 
(1.1) [8] suggests it is more natural to think of a cancellation between n|j,!f and twice the self-energy (3.13), 
and between Hfj,^ and twice the vertex correction (3.16). For the analogous diagrams in the Gross-Neveu 



case (note the constant Ad is different): 

nMr^^ - _„(fc^9SW(fc) _ (k^)i-Ud~2) Cd A. 
^'~ Ad dm ~ Ad d Nf^'k' 

Ad i-g/^) Ad Nf k 

In fact, what has occured is a cancellation between the subdivergences of figure 4 and the diagrams which 
would result from inserting local counterterms arising from the divergences of eqns. (3.13) and (3.16) in the 
leading order vacuum polarization figure 1 - factors of two come because there are two fermion lines and 
two vertices to correct. 
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In a treatment of renormalization which proceeds by subtraction of divergent parts in an ordered fashion 
- the BPHZ scheme - this cancehation is transparent. In the presentation here I have chosen a physically 
more intuitive approach - rescaling the fields and coupling parameters in the bare Lagrangian to keep Green 
functions finite at each order - but of course the two schemes are completely equivalent. However, it is 
important to note that the argument depends on a novel application of Weinberg's theorem [20]. The 
theorem states, in effect, that a IPI graph such as IlJ^J, with positive degree of divergence lo will be balanced 
by an overall countcrterm which is polynomial in external momentum (in this case 0(A''^^fc*')), provided 
that its subdivergences are first subtracted, in this case by applying the appropriate vertex and self-energy 
corrections. These subdivergences may be non-polynomial in momentum, but the theorem guarantees that 
the cancellation will be exact. The crucial point is that Weinberg proved the theorem for a wide class of 
integrands; the usual integrals built from the standard Feynman propagators proportional to {p^ + m?)^^, 
(i'^ + m)^^ etc. form just one class. There is no general requirement that the propagators be algebraic in 
the internal momenta, and integrals including the resummed propagator ex k'^^'^ arc also included. ^ 

To the reader familiar with the calculation of the two-loop vacuum polarization in QED4 (eg [21] ch. 
8), it is worth making a further comment. The analogue of the k'^^^ In A diveregnces are not (as one might 
first think) the fc^ In A terms which lead to a physical charge renormalization - these are polynomial in k 
and "belong" to the diagram as a whole - but rather the non-polynomial InfclnA divergences, which in 
general are cancelled by countcrterm subtractions, and in the particular case of QED cancel between the two 
diagrams (4a) and (4b). The discontinuous behaviour as d ^ 4 illustrates why this is a critical dimension 
for the model. 

So, we see that in a sense the main achievement of this paper is simply the verification of a peculiar case 
of Weinberg's theorem. There appears to be no fundamental obstacle to formulating a proof of rcnormaliz- 
ability of the 1/iV/ expansion for four-fermi theories to all orders. The only complication arises, as we have 
seen, because graphs with different numbers of loops arise at a given order, which means that countcrterm 
subtractions oi the same order must be applied to yield a finite result, or in the language of [8], non-trivial 
cancellations between divergent graphs at the same order must occur. As we argued in the last section, this 



^ In fact, one requirement of Weinberg's original proof is that the momenta are defined in an integer- 
dimensioned vector space; the extension to non-integer d has not been established to my knowledge. There 
are two responses: either set d = 3 at this stage to yield a "physical" theory, or note that it is always possible 
to route a simple loop momentum through any internal auxiliary line, in which case it may be possible to 
analytically continue the integrand to d = 3 where the theorem holds. We shall not pursue these rather 
abstract issues further. 
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does not occur in the l/Nf expansion of the nonhnear signia niodcL 

We deduce that if the expansion is rcnormahzable, the logarithmic corrections to fermion and vector 
propagators and the vertex always cancel at each order in l/Nf. An important physical consequence, which 
also follows from Weinberg's theorem, is that in both Gross-Nevcu and Thirring models, the amplitude 
for four fermion scattering receives no l/Nf corrections in the deep Euclidean limit [8]. In both cases it 
assumes a universal form Ad/{Nfk'^^'^), which thus characterizes an interacting ultraviolet fixed point of the 
renormalization group. Both models resemble each other at high energies. At low energies they differ, and 
now we return our focus to the Thirring model. In section II we saw that at leading order the coupling g is 
completely unconstrained, and that the model can be formulated cither as weakly or strongly coupled. After 
radiative corrections, this may no longer be true. Due to the mass renormalization (3.14), the fermion mass 
operator acquires an anomalous dimension of 0{1/Nf). The result is that 

M °^ (to j ^ (4-2) 

that is, for fixed physical mass to,, the bare mass M must be tuned to zero as the cutoff is removed. 
Accordingly the ratio g^ /A'P^'^ , which governs whether the model is weakly or strongly coupled at leading 
order, must grow small. However, since the low energy nature of the model is characterized by the ratio 
My /to,, it will be necessary to compute 0{1/Nf) corrections to My to determine whether the model is driven 
to strong or weak coupling at higher orders. It would also be interesting to test the stability of (4.2) under 
corrections of 0{\/N'i). 

Finally we speculate on the relevance of this model to strongly-coupled QED, both in 3 and 4 dimensions. 
This paper has been concerned exclusively with l/Nf perturbation theory. In [19,23], the leading order vector 
propagator (2.15) was used in the Schwinger-Dyson equation to solve for dynamical fermion mass generation 
self-consistcntly. The result, oc exp(— A^^), is non-perturbative in l/Nf. It is suggested that the solution may 
shed light on spontaneous chiral symmetry breaking in QEDg , which is suspected to be critically dependent 
on Nf [24]. Since QEDg is super-renormalizable, the continuum limit is thought to exist in the limit of weak 
coupling; so far we can reach no conclusion for the Thirring model. It is also worth noting that since all 
0{1/Nf) corrections to the vector propagator are UV finite, then in the asymptotic regime the quenched 
approximation (which must be made to solve the Schwinger-Dyson equation) is exact. This is therefore a 
model with many similarities to QED in which charge screening is naturally switched off at short distences 
- an effect which must be postulated in studies of a non-trivial fixed point of QED4 due to fermion mass 
generation at strong coupling. This deserves further study. Another possibility is to introduce both scalar 
and vector current interactions with independent couplings, to generate a fermionic analogue of the gauged 
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Nambu - Jona-Lasinio model, which would be renormalizable in rf G (2,4). It would then be interesting to 
examine the limit d ^ 4_ and compare the results with other approaches [25]. 
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Figure 1 

Leading order contribution to the vector auxiliary two-point function 
Full lines represent fermions, wavy lines the vector auxiliary 




Figure 2 

0{\/Nf) contribution to the fermion self-energy 
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Figure 3 

0{1/Nf) contribution to the vertex 
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Figure 4a 

0(1/Nf) contribution to the vector two-point function 



1^ 




Figure 4b 

0{1/Nf) contribution to the vector two-point function 
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Figure 5a 

Diagram representing J^i^api^) 
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Figure 5b 

Diagram representing J^^^ai^) 
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